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Abstract. We re-analyzed soft X-ray data collected with the ROSAT and ASCA observatories on a candidate neutron 
. star (NS) near the center of the supernova remnant PKS 1209-52. We fitted the observed spectra with NS atmosphere 
models. The hydrogen atmosphere fits yield more realistic parameters of the NS and the intervening hydrogen column 
. than the traditional blackbody fit. In particular, for a NS of mass 1.4 Mq and radius 10 km, we obtained the NS surface 
■ temperature T^s = (1.4 - 1.9) x 10^ K and distance d = 1.6 - 3.3 kpc versus T = (4.2 - 4.6) x 10^ K and (implausible) 
d — 11 — 13 kpc for the blackbody fit, at a 90% confidence level. Our fits suggest that the surface magnetic field is either 
very weak, B < 10^" G, or it exceeds ~ 2 x 10^^ G. The hydrogen column density inferred from the atmosphere fits, 
T— I ' = (0.7 — 2.2) X 10^^ cm^^, agrees fairly well with independent estimates obtained from UV observations of nearby 
^ stars, radio data, and X-ray spectrum of the shell of the supernova remnant, whereas the blackbody and power-law 
l> , fits give considerably lower and greater values, uh = (0.2 — 0.4) x 10^^ and (5.2 — 7.0) x 10^^ cm ^, respectively. The 
' inferred NS surface temperature is consistent with standard NS cooling models. 
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2 1- Introduction 



In addition to about 15 radio pulsars associated with supernova remnants (Gaensler & Johnston 1995), several radio 
silent isolated neutron star (NS) candidates within SNRs have been observed with the X-ray observatories HEAO A, 
Einstein, EXOSAT, ROSAT, and ASCA (see Caraveo, Bignami & Triimper 1996 for a review). These objects have not 
been detected outside the X-ray range, and their X-spectra resemble blackbody (BB) spectra with temperatures of a 
5— i ' few million kelvins. If they are indeed thermally emitting NSs, the analysis of their radiation provides an opportunity 
to study thermal evolution of NSs of ages 10^ — 10^ yr, which is important for elucidating the properties of the 
superdence matter in NS interiors. 

One of the most convincing examples of such objects is the point-like source IE 1207.4-5209 within the barrel- 
shaped radio. X-ray, and optical SNR PKS 1209-52 (also known as G269. 5-1-10.0). From the analysis of radio and 
optical observations of this SNR, Roger et al. (1988) estimated its age ~ 7000 yr, with an uncertainty of a factor 
of 3, and concluded that the remnant is in an adiabatic expansion phase. The distance to PKS 1209-52 is not well 
determined — estimates in the range 1.1 — 3.9 kpc were suggested (Milne 1979; Mills 1983). Estimates of the interstellar 
hydrogen column density from the radio and optical data yield n/f 21 = n^/iW^^ cm~^) ^ 1.0 — 1.8 (see Roger et 
al. 1983 and Kellet et al. 1987 for references), consistent with a distance d ^ 1 — 2 kpc. 

After the first X-ray detection of PKS 1209-52 with HEAO A-1 (Tuohy et al. 1979), the point source IE 1207.4- 
5209 was discovered with the Einstein observatory (Helfand & Becker 1984), 6' off-center the 81' diameter SNR. Matsui, 
Long & Tuohy (1988) concluded that its spectrum can be interpreted as a BB spectrum with an apparent temperature 
T = 1.4 X 10^ K (assuming n/f 21 = 3.2, as obtained from the HEAO A-1 observations). The Einstein High Resolution 
Imager (HRI) observations showed a lack of diffuse X-ray emission around IE 1207.4-5209 (later confirmed by the 
ROSAT HRI observations), which greatly simplifies the analysis of the point source radiation. From the analysis of 
EXOSAT observations, Kellett et al. (1987) estimated the BB temperature of the central object T = 1.8 x 10^ K (at 



2 



Zavlin et al.: Neutron Star in PKS 1209-52 



nH,2i = 1-1), which impUes an emitting area of a radius (3 — 4) d2 km, where d2 = d/ (2 kpc). Applying the Raymond 

& Smith (1977) linc-cmission model, Kellet et al. estimated also an average SNR temperature, ^ 1.7 x 10^ K, and the 
hydrogen column towards the SNR, nH,2i — 1-4, consistent with the values obtained from radio and optical data. 

Observations of PKS 1209-52 and its central source with ROSAT and ASCA have further supported the NS 
hypothesis for IE 1207.4-5209. Mereghctti, Bignami & Caraveo (1996; hereafter MBC96) showed that the ROSAT data 
on IE 1207.4-5209 can be interpreted as blackbody emission of T ~ 3 x 10^ K from an area with radius R ~ 1.5 d2 km. 
The authors noticed that this temperature is too high to be explained in the framework of a cooling NS with age of 
PKS 1209 52 (~ 10^ yr). The hydrogen column inferred from the BB fit, nH.21 ^ 0.4, is 3 4 times lower than the 
estimate given by Kellett et al. (1987). From observations at 4.8 GHz, MBC96 found an upper limit of ~ 0.1 mJy on 
the radio flux from IE 1207.4-5209. They also set a deep limit of F > 25 for an optical counterpart in the Einstein 
HRI error circle, that supports the hypothesis that IE 1207.4-5209 is indeed an isolated NS. 

Vasisht et al. (1997; hereafter Val97) have recently analysed ASCA observations of IE 1207.4-5209. They found 
that each of the three spectra obtained with the ASCA detectors can be fitted by a BB spectrum consistent with that 
obtained from the analysis of the ROSAT data. The hydrogen column is poorly restricted in the ASCA spectral fits due 
to low detector sensitivities at photon energies below 0.5 keV. Using a fit with a Raymond-Smith model at fixed cosmic 
abundances to the ROSAT data on the SNR shell, Val97 estimated the remnant temperature ~ (1.9 — 2.3) x 10^ K and 
the foreground column density nH.21 ^ 0.6 — 0.9. They attributed the difference in uh between the NS and SNR fits to 
either separate lines of sight or the large column density that the shell X-rays can encounter in the SNR postshocked 
gas. 

Based on the results of the blackbody analysis, both MBC96 and Val97 concluded that IE 1207.4-5209 is an 
isolated NS with hot spots on its surface, aligned with the magnetic poles. Val97 suggested that the spots are heated 
either by dissipative heating in the NS interior or by the bombardment of polar caps by relativistic particles from the 
NS magnetosphere if IE 1207.4-5209 is an active pulsar. However, the former hypothesis can hardly explain the small 
sizes of the hot spots, even with allowance for large anisotropy of thermal conductivity of the magnetized NS crust. 
The latter heating mechanism is also in doubt because of the absence of radio and 7-ray emission from IE 1207.4-5209. 
Moreover, both the ROSAT and ASCA data did not reveal pulsations, which may not be consistent with the presence 
of the hot spots unless the magnetic and rotation axes are coaligned. 

Although it looks very plausible that IE 1207.4-5209 is a thermally emitting isolated NS, the BB interpretation 
adopted by previous authors leaves several controvertible points. On the other hand, thermal radiation emitted by NS 
atmospheres may significantly differ from the BB radiation. Therefore, to resolve the inconsistencies following from 
the simplified BB interpretation, it is natural to employ more realistic models of NS radiation. Here we present a 
combined analysis of the ROSAT and ASCA data based on NS atmosphere models (Pavlov et al. 1995, and references 
therein). These models have been applied successfully to the interpretation of the soft X-ray radiation from, e.g., 
the Vela pulsar (Page, Shibanov & Zavlin 1996) and the brightest millisecond pulsar J0437-4715 (Zavlin & Pavlov 
1997; Pavlov & Zavlin 1997). These examples show that fitting the soft X-ray pulsar spectra with hydrogen or helium 
atmosphere models always results in lower effective temperatures and greater emitting areas (or smaller distances) 
than those obtained from the BB fits. 

We show in §2 that, indeed, the model atmosphere fits of the X-ray radiation from IE 1207.4-5209 yield an NS 
surface temperature compatible with NS cooling models, and they do not require hot spots on the NS surface. Moreover, 
the hydrogen column density inferred from this interpretation is in excellent agreement with that obtained from our 
fits of the SNR X-ray radiation as well as with independent estimates of nn for stars in the vicinity of IE 1207.4-5209 
(§3). These results warrant application of the same approach to other similar objects and enable one to obtain reliable 
estimates of surface temperatures of NSs of different ages. 

2. Observations and results 

We used archival ROSAT and ASCA data on IE 1207.4-5209 (see MBC96 and Val97 for a detailed description). 
Four ROSAT Position Sensitive Proportional Counter (PSPC) pointings (obtained in 1993 July, with exposure times 
« 5 — 6 ks each — see Table 2 in MBC96) were centered at different positions (off-axis angles between 16' and 33'). We 
checked that the spectra of IE 1207.4-5209 extracted from each of the four data sets and properly corrected for the 
PSPC response are consistent with each other. For the combined analysis with the ASCA data we chose the ROSAT 
PSPC observation with the minimum off-axis angle 16' (exposure time 4.97 ks). We extracted a raw spectrum from a 
region with a radius r ~ 75" centered on the source. The (off-axis) source count rate is 0.145 ± 0.005 s"^. We binned 
this spectrum into 24 spectral bins in the 0.1-2.4 keV energy range for further analysis. The ROSAT HRI counts 
for each of the three on-axis pointings (obtained in 1992 August and 1994 July; total exposure time « 14.5 ks) were 
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selected from the circles of r ~ 22". The corresponding source count rate is 0.064 ± 0.003 s ^, in agreement with 
MBC96. 

The ASCA data were obtained in 1994 July, with two Solid State Imaging Spectrometers, SISO and SISl, in the 
bright mode, and two Gas Scintillation Imaging Spectrometers, GIS2 and GIS3, in the pulse-height mode. The SISl 
image of IE 1207.4-5209 is too close to the edge of the CCD chip for a reliable spectral analysis. The SISO source 
counts collected in w 21 ks (with high and medium bit rates) were extracted from the circle of r ~ 4', the total source 
count rate is 0.0556 ±0.0017 s~^. The SISO spectrum was binned into 38 spectral bins in the 0.5 — 10.0 keV range. The 
GIS2 and GIS3 counts were extracted from circles of r ~ 9' (exposure w 19 ks for the high plus medium bit rates). 
The source count rates are 0.0425 ± 0.0018 (GIS2) and 0.0533 ± 0.0025 s'^ (GIS3). The spectra were binned into 40 
and 45 spectral bins in the 0.3 — 5.0 keV range for GIS2 and GIS3, respectively. The total source count rates for all 
the three ASCA instruments are somewhat lower than those reported by Val97, perhaps because of different screening 
criteria we applied. 

As a first step of our data analysis, we fitted the count rate spectra simultaneously for the four instruments, 
ROSAT PSPC and ASCA SISO, GIS2, and GIS3, with the traditional power-law and BB models. Figure 1 shows 
the results of these fittings. The power-law fit yields a photon index 7 = 5.2(-|-0.4, —0.3) and a hydrogen column 
density, nH,2i = 6.0(+1.0, —0.8) (minimum xt = 1-47, v = 144; the uncertainties hereafter are at a 90% confidence 
level). These values of 7 and uh are better constrained than those obtained by MBC96 from the PSPC data alone, 
7 = 3.9(-|-1.5, — 1.4) and nH,2i = 3.0(-|-2.5, — 1.8) (minimum xt = 1-61, = 10). The photon index is unusually 
large in comparison with typical values, 7 ~ 1.5 — 2.5, observed from, e.g.. X-ray emitting radio pulsars (Becker & 
Triimper 1997), and the hydrogen column density too much exceeds the values obtained by independent measurements 
(see below). These discrepancies, together with the too high minimum xt, imply that the power-law interpretation is 
inadequate. 

The BB temperature (as measured by a distant observer). T°° = 3.12(+0.15, —0.11) x 10^ K, the apparent radius 
of the emitting area, Ra — 1.21(-|-0.10, -0.13) ^2 km, and the bolometric luminosity, Lg^i = 0.93(-|-0.19, -0.11) x 
10^^ (^2 ergs~^, obtained from the combined fitting (Fig. 1) are very close to the values from the separate fits of 
the ROSAT and ASCA data by MBC96 and Val97. The hydrogen column density in the combined fitting, nH.21 = 
0.28(-|-0.15, — 0.12), is compatible with (albeit slightly lower than) the value of nH,2i = 0.3 — 0.5 inferred from the 
separate ROSAT PSPC spectrum. Although the quality of the BB fit is better than of the power law fit, the uh 
values look surprisingly low in comparison with (1.0 — 1.8) x 10^^ cm~^ obtained from independent estimates (see 
Introduction) . 

Since thermal NS spectra, as well as spectra of ordinary stars, cannot exactly coincide with the BB spectra because 
of the effects of radiative transfer in the emitting layers, it is natural to compare them with more realistic model 
spectra of NS atmospheres (e.g., Pavlov et al. 1995). Strongest deviations from the BB spectrum are expected if 
the NS surface is covered with a hydrogen atmosphere. Hydrogen may appear at the NS surface as a result of, e.g., 
accretion of interstellar matter or post-supernova accretion of a fraction of the ejected envelope. Due to the huge 
surface gravity (typical gravitational acceleration ~ 10^^ — 10^^ cm s^'^), heavier chemical elements sink down in 
deeper layers and do not affect properties of emitted radiation, whereas hydrogen remains at the surface. The shape 
of the spectrum emitted by an atmosphere depends on the strength B of the NS surface magnetic field. The lack of 
pulsations does not allow one to estimate B using period and its derivative, as is usually done for pulsars. The absence 
of statistically reliable features in the observed spectra, which could be associated with an electron cyclotron line at 
Ebc = 11.6(i3/10"'^^ G) keV, prevents one to make a direct estimation of B. Therefore, one has to try the atmosphere 
models with both low {B < 10^" G) and high magnetic fields. (In the former case, the field does not affect properties 
of emergent radiation [Zavlin, Pavlov & Shibanov 1996], so that we can merely put B = 0). In the case of strong 
nonuniform field {B > 10^^ G), the NS should have a nonuniform temperature distribution along its surface because of 
the high anisotropic thermal conductivity in the NS crust (e.g., Shibanov & Yakovlev 1996). However, since we have 
no information about the field geometry, we assume that the magnetic field has the same strength and is directed 
radially everywhere at the surface, and the effective temperature is uniform. This assumption reduces the number of 
fitting parameters and can be considered as a reasonable first approximation for investigating the atmosphere effects. 

Figure 2 shows results of fitting of the combined ROSAT and ASCA data with hydrogen atmosphere models 
for three fixed values of the magnetic field. In these fits we assume canonical values for the NS mass and radius, 
M = I.4M0 and i? = 10 km, and consider the distance as a fitting parameter. The models with B = (left panels) 
and 10^^ G (right panels) result in close values of the distance, effective temperature at the NS surface and bolometric 
luminosity, d = 2.3(+1.0, -0.7) kpc. Toff = 1.63(-|-0.27, -0.19) x 10^ K, iboi = 5.0(+4.3, -1.8) x lO^^ ergs-\ and 
d = 2.3(-M.0,-0.8) kpc, Teff = 1.65(-h0.20, -0.17) x 10^ K, Lboi = 5.3(-h3.0, -2.0) x 10^3 ergs'^, for S = and 
10^3 G, respectively. Since the nonmagnetic spectra are softer at lower energies, the hydrogen column density at 
B = 0, nH,2i = 1.0(+0.5, —0.6), is lower than for the strong magnetic field, nH,2\ = 1.5(+0.7, —0.6). We checked that 
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folding the models within the 90% confidence region with the ROSAT HRI response yields count rates compatible 

with those observed. When the field strength varies in the range 5 x 10"'^^ ^ B <5 x 10^^ G, the fitting parameters 
remain approximately the same because the model spectra are almost insensitive to the B value in the corresponding 
domain of energies and effective temperatures. When B exceeds ~ 5 x 10^^ G, the proton cyclotron line centered at 
Ebp = 0.063(-B/10"'^^ G) keV (cf. Bezchastnov et al. 1996) enters the SIS energy range, which makes the fits statistically 
unacceptable. (This line moves above ~ 5 keV, a maximum energy where the NS fiux is still above the background, 
at superstrong magnetic fields, B > 8 x 10^^ G, for which the models we used here are not directly applicable.) When 
B falls below ~ 5 x 10^^ G, the low-energy wing of the electron cyclotron line gets into the ASCA range, and the 
model spectra become softer at E >2 keV. As a result, the confidence contours in the nn-d and nH-T^s planes move 
to greater d and Teg, and lower ur, towards the BB contours. An example is shown in the middle panels of Fig. 2 for 
B = 2 X 10^^ G, for which the best-fit distance is about 50% larger than at i? = lO""^"^ G. When the field is lower than 
~ 5 X 10^^ G, but greater than ~ 1 x 10^° G, the atmosphere model fits become statistically unacceptable because the 
core of the electron cyclotron line gets into the ASCA/ ROSAT range. 

The atmosphere models depend not only on _B, but also on the NS mass and radius which determine the gravitational 
acceleration (one of our model parameters) and the gravitational redshift factor = (1 — 0.295M«/i?io)^/^, where 

= M/Mq, i?io = R/10 km. To illustrate this effect, we present in Fig. 3 the best-fit parameters at B = 10^^ G in a 
wide range of R and M allowed by equations of state of the NS matter. Although the effective temperature at the NS 
surface. Teg, varies by about ±20% in the allowed R-M domain, the apparent effective temperature (as measured by 
a distant observer), = grT^a, remains almost constant, (1.23— 1.31) x 10^ K, because the redshift is compensated 
by the change of the unredshifted NS spectrum (it softens in the Wien tail with increasing M and decreasing R at 
given Teff). Owing to the approximate constancy of T^, the apparent and "true" bolometric luminosities depend on 
R and M as ig^j oc R^ = R^ /g^ and Lboi = P^^-^boi ^ ^ l9r- The latter dependence explains the non-monotonous 
behaviour of L\^oi at higher M. The best-fit d and uh are almost independent of M at higher i?, when is close to 
1. At lower R the variations are stronger, albeit within the statistical uncertainties of these parameters (cf. Fig. 2). 
The distance inferred at assumed R and M can be approximately described, in the allowed mass-radius domain, by a 
linear equation: d = 0.5M* — 0.4 -I- (2.5 — 0.3M*)i?io (in kpc); it grows with R faster than for the BB interpretation. 
Notice that if d is determined more accurately in future observations of the SNR and its central source, this equation 
would delimit a band in the M-R plane constraining equation of state of the NS matter. 

Since the hydrogen atmosphere fits yield uh greater than the BB fit by a factor of 3-5, it is important to estimate 
this parameter independently. For this purpose, we fitted the SNR emission with various models for thermal plasma 
radiation. We extracted the SNR spectrum from a bright region of the ROSAT PSPC image of the remnant shell within 
the r ~ 5' circle centered at a(2ooo) = 12'^11™56'?0, (5(2ooo) = -52°36'36". Three models available with the XSPEC 
package, vraymond, vmeka, and vmekal (Raymond-Smith, Mewe-Gronenschild-Kaastra, and Mewe-Kaastra-Liedahl 
models with variable abundances) give satisfactory fits (x^ < 1-0) with the column density to the SNR in the range of 
(1.1 — 1.9) x 10^"'^ cm~^, consistent with the hydrogen atmosphere fits of IE 1207.4-5209, but certainly in excess of the 
BB fits. The inferred plasma temperature is (1.9 — 2.1) x 10^ K. These results were obtained with a moderate excess 
of abundances of Al and Si whose emission lines are prominent in the SIS spectra of the SNR shell. Abundances of 
other elements are close to standard values. The SNR parameters inferred from our fittings are consistent with those 
obtained by Kellett et al. (1987) from the EXOSAT data. The lower hydrogen density, nii,i\ = 0.6 — 0.9, obtained by 
Val97, is likely associated with fixed cosmic abundances adopted, that resulted in lower statistical quality of their fits 
ixl = 1-8). 

3. Discussion and conclusions 

We have shown that the X-ray spectra of the isolated NS IE 1207.4-5209 in the center of the SNR PKS 1209-52 
observed with the ROSAT PSPC and HRI and ASCA SISO, GIS2 and GIS3 can be interpreted as thermal radiation from 
the hydrogen-covered, uniformly heated surface of the NS. The proposed interpretation resolves all the inconsistencies 
which follow from the BB interpretation of radiation of this object. 

First, the range of distances, d = 1.5 — 3.2 kpc, inferred for the standard NS mass M = I.4M0 and radius 
i? = 10 km, is within the (conservative) limits, 1.1 — 3.9 kpc, obtained from radio observations of the SNR. This 
range shifts to lower distances with decreasing R and M: e.g., d = 1.2 — 2.6 kpc for M = l.OM©, R = 8 km. Thus, 
the assumption about hot spots on the NS surface, which is difficult to reconcile with the lack of pulsar activity in 
IE 1207.4-5209, is superfiuous in our model. 

Second, the effective temperature Tes = (1.4 - 1.9) x 10*^ K {T^ = (1.0 - 1.6) x 10*^ K), obtained from the fittings 
with our hydrogen atmosphere models, matches well to a number of the NS cooling models (see Fig. 4 and discussion 
below). On the contrary, it is difficult to reconcile the results of the BB interpretation with models of the NS cooling. 



Zavlin et al.: Neutron Star in PKS 1209-52 



5 



Third, the inferred hydrogen column density towards the IE 1207.4-5209, nH,2i = 0.9 — 2.2 and 0.5 — 1.5 for 

high and low surface magnetic fields, respectively, agrees fairly well with 71^,21 = 1-1 — 1-9 obtained from fitting the 
spectra of the SNR shell. These estimates are consistent with the nn values obtained by Kellett et al. (1987) from 
the EXOSAT observations of PKS 1209-52. Moreover, measurements of uh in UV observations of several stars in the 
vicinity of IE 1207.4-5209 (/ = 296.5°, b = 9.9°) yield the hydrogen columns in virtually the same range (Fruscionc et 
al. 1994): ni/,21 = 1.1 - 1.5 for HD112244 {I = 303.6°, b = 6.0°, d = 1.85 kps), ni/,21 = 1.4 for HD115842 (/ = 307.1°, 
b = 6.8°, d = 1.87 kps) and HD111822 {I = 303.1°, b = 10.2°, d = 2.53kps).'On the contrary, the BB fit gives 
n'H,2i = 0.15 — 0.45, clearly incompatible with all the other independent estimates. 

The fact that IE 1207.4-5209 shows no pulsar activity can be explained, in addition to the trivial explanation of 
unfavorable orientation of its magnetic and rotational axes, by slow rotation or/and low magnetic field of the NS, so that 
it falls below the "death line" on the P-P diagram, P/P^ = 2 x 10^^^ . For the simplest model of the magnetic 
dipole radiator, the pulsar period increases with time t as P = {P^ + 2^o^)^/^ where to = {Stt'^ /3){B^R^ /Ic^) = 
0.97 X 10-13^2^^6^ j-i g J ^ lO'^^hr, g cm^ is the moment of inertia of the NS and B13 = B/{10^^ G). The condition 

that the NS is below the death line can be written as io/P* < 2 x 10"^'^ s'^, or P > SABl^^Rld^I^^^^^ s. Hence, if 
B is in the range of strong magnetic fields allowed by our fitting, 0.2 < B13 < 8, we may expect that the NS rotates 
slowly, P > 1 — 2 s. If faster pulsations are discovered in future X-ray observations of IE 1207.4-5209, it would indicate 
that B is in the other range compatible with our fits: B < 10^° G. It is worth noting that if the NS was born as a 
pulsar at an appreciable distance from the death line, only an enormous mag netic field, B > 5 x 10^'^ G at standard NS 
parameters, could decelerate its rotation so that the pulsar would "die" at the relatively young age, t ~ 10^ yr. Hence, 
discovery of slow pulsations would mean that the NS was born slowly rotating unless its magnetic field is superstrong. 

Since we have estimated the effective temperature of the NS, it is interesting to compare it with what is predicted 
by various cooling models. Figure 4 shows examples of NS cooling curves from Van Riper, Link & Epstein (1995) for 
three equations of state (EOS), stiff (Pandharipande-Smith [PS]), intermediate (Priedman-Pandharipande [FP]) and 
soft (Baym-Pethick-Sutherland [BPS]), and two interior compositions resulting in slow and fast cooling (solid and 
dashed lines, respectively). Cooling curves without additional heating are depicted by thick lines, whereas thin lines 
show cooling curves for two models, proposed by Epstein & Baym (EB) and Alpar, Cheng & Pines (ACP), for pinning 
of the superfiuid vortices to the crust lattice. They correspond to strong and weak frictional heating associated with 
the dissipation of energy of differential rotation between the NS crust and superfiuid interior (only the EB model is 
available for the BPS EOS). In the same picture we plot boxes corresponding to the inferred atmosphere and BB 
temperatures of IE 1207.4-5209, assuming its age in the range (0.3 — 2) x 10"* yr, and similar boxes for the Vela pulsar 
whose X-ray radiation was investigated in terms of the hydrogen atmosphere models by Page et al. (1996). (We adopted 
(1 — 4) x 10* yr for the Vela age, between the conventional characteristic age P/2P and an upper limit estimated by 
Lyne et al. 1996). We see that the BB temperature of IE 1207.4 5209 is well above the values predicted by all these 
cooling models, whereas the BB temperature of the Vela pulsar is compatible only with the standard (slow) cooling 
model supplemented by strong heating for the stiff EOS. The effective temperature of IE 1207.4-5209 obtained with 
the atmosphere fits is compatible, given the poorly known age, with the slow cooling models for all the three EOS, at 
moderate or no heating. For the stiff EOS, it is also compatible with the fast (quark) cooling accompanied by strong 
frictional heating. However, the large NS radius for this EOS would mean a distance ~ 2.5 — 5.2 kpc (cf. Fig. 3), 
uncomfortably large in comparison with conventional estimates. The same cooling curve (fast cooling with the EB 
heating for the stiff EOS) is the only one which goes through both the IE 1207.4-5209 and Vela pulsar boxes obtained 
from the atmosphere fits. This, however, does not necessarily mean that slow cooling models or other EOS or other 
heating rates are excluded. First, the ROSAT PSPC spectrum of the Vela pulsar, used in both BB and atmosphere 
fits, was not observed directly because the pulsar was not resolved from a surrounding mini-nebula of 2' diameter, 
which makes the results of the spectral fits less certain. Second, one cannot exclude, in principle, that the NS of the 
Vela pulsar differs from IE 1207.4-5209 (e.g., the mass of the former may be greater, that could lead to an "exotic" 
interior composition associated with an enhanced neutrino luminosity and accelerated cooling). Finally, there exist 
many more cooling models than shown in Fig. 4, and some of them may satisfy both the IE 1207.4-5209 and Vela 
pulsar temperatures. For instance, a strong neutrino emission induced by the nucleon Cooper pair formation process 
(see Page 1997, and references therein), which was neglected in most of the previous NS cooling models, may result 
in great variety of cooling curves, depending on the (unknown) parameters of the nucleon pairing (Yakovlev 1997). 
For a more detailed comparison of the inferred temperatures with the cooling models, it would be very important to 
evaluate more accurately the distance to PKS 1209-52 (which would constrain the NS radius and hence the EOS) and 
the age of this object. 

There exist a number of other radio silent isolated NS candidates whose observational manifestations are similar 
to those of IE 1207.4-5209. The most convincing example is lE/RXJ 0820-4247 in the SNR Puppis A, whose BB 



6 



Zavlin et al.: Neutron Star in PKS 1209-52 



temperature, « 3 x 10^ K, and radius, w 2 km (Petre, Becker & Winkler 1996), virtually coincide with those of 
IE 1207.4-5209. Another example with similar properties is RXJ 0002+6246 in the SNR G117.7+0.6 (Hailey & Craig 
1995), whose identification is, however, less certain because of its faintness (likely, due to a larger distance and greater 
hh)- We expect that applying the NS atmosphere models to the analysis of such objects will allow us to evaluate their 
radii and effective temperatures, and to constrain their magnetic fields. 

One more radio silent NS candidate, IE 1613-5055 in the center of the SNR RCW 103 (Tuohy & Garmire 1980), 
appears as a different kind of object — the BB fitting of its spectrum yields a considerably higher temperature, 
7 X 10^ K, at a few times smaller radius of the emitting region (Gotthelf et al. 1997). However, this object is deeply 
immersed in the remnant diffuse emission so that it is hard to separate its spectrum from that of the SNR, with the 
limitied spatial resolution of ASCA. We expect that the forthcoming AXAF and XMM missions would resolve the 
point source and provide its spectrum suitable for the detailed analysis. 
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Figures 

Fig. 1 68%, 90% and 99% confidence contours for the power-law and blackbody fittings to the ROSAT PSPC and ASCA 
SISO, GIS2 and GIS3 spectra of IE 1207.4-5209. The fines in the right panel correspond to constant values of radius of emitting 
area at a distance d = 2 kpc. 

Fig. 2 68%, 90% and 99% confidence contours for the NS hydrogen atmosphere fittings for M = IAMq, R = 10 km, and 
three values of the surface magnetic field B to the combined ROSAT and ASCA data. 

Fig. 3 Dependences of the best-fit parameters for the atmosphere fittings on NS radius R at different values of the NS mass, 
M/Mq — 0.8, 1.0, . . . 2.0. The thick dashed lines delimit the ranges of the apparent temperature and luminosity (upper and 
lower dashed curves correspond to M/Mq = 2.0 and 0.8, respectively). 

Fig. 4 Comparison of the surface temperatures inferred from the BB (empty boxes) and hydrogen atmosphere (hatched 
boxes) fittings of the spectra of IE 1207.4-5209 and the Vela pulsar with NS cooling curves for three EOS of the NS matter, 
two NS compositions, and different models for additional frictional heating. The cooling curves are taken from Van Riper et 
al. (1995), where more detailed explanations can be found. 
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